function M = elipseMask(SizeX, SizeY, el )
%ELIPSEMASK vrati elipsu (plnou) jako binarni obrazek
%           SizeX, SizeY ... rozmery vysledneho obrazku
%           el - elipsa ve forme:
%              el(1)x^2 + el(2)xy + el(3)y^2 + el(4)x + el(5)y + el(6) = 0


if el(1) < 0 | el(3) < 0
    el = -el;
end
[X, Y] = meshgrid(1:SizeX, 1:SizeY);
M = (X.*(el(1)*X+el(2)*Y+el(4))+Y.*(el(3).*Y+el(5))+el(6)) <= 0;

end

